oe- /S9, z f/ 


NASA Contractor Report 159294 


NAS A-CR-l 59294 

Hgo GO / "SDlQI 


VISCOELASTIC ANALYSIS OF ADHESIVELY 
BONDED JOINTS 


F. Del ale and F. Erdogan 


LEHIGH UNIVERSITY 
Bethlehem, Pennsylvania 18015 


NASA Grant NGR 39-007-011 
March 1980 


ivinsA 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton. Virginia 23665 




NF01094 


VISCOELASTIC ANALYSIS OF ADHESIVELY 
BONDED JOINTS* 

by 
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Abstract 

In this paper an adhesively bonded lap joint is analyzed by 
assuming that the adherends are elastic and the adhesive is linearly 
viscoelastic. After formulating the general problem a specific example 
for two identical adherends bonded through a three parameter visco- 
elastic solid adhesive is considered. The standard Laplace transform 
technique is used to solve the problem. The stress distribution in 
the adhesive layer is calculated for three different external loads 
namely, membrane loading, bending, and transverse shear loading. The 
results Indicate that the peak value of the normal stress in the 
adhesive is not only consistently higher than the corresponding shear 
stress but also decays slower. 


(*) This work was supported by NASA-Langley under the Grant NGR 39“007~011 
and by NSF under the Grant ENG 78-09737* 
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1. Introduction 


In Its simplest form an adhesively bonded structure consists of 
three, components of different mechanical properties, namely the adhe- 
sive 'and the two adherends. Because of the nonhomogeneous nature and 
of the geometrical complexity of the medium, even for the linearly 
elastic materials the exact analytical treatment of the problem regard- 
ing the stress analysis of the structure Is, in general, hopelessly 
complicated. The existing analytical studies are, therefore, based 
on certain simplifying assumptions with regard to the modeling of the 
adhesive and the adherends. The adherends are usually modeled as an 
Isotropic or orthotropic membrane (e.g., [l]), plate (e.g., [2,3]), or 
elastic continuum (e.g., [A, 5]). The primary physical consideration • 
used in the selection of a particular model is generally the ratio of 
the thickness of the adherend to the lateral dimensions of bond region. 
For example, for adherends with a very small relative thickness the 
bending stiffness may be neglected whereas if the thickness of the 
adherend is not small even the plate assumption may be erroneous. As 
for the adhesives, generally the thickness variation of the stresses 
is neglected and the adhesive layer is modeled as a linear shear or a 
tension-shear spring. 

In most applications of structural adhesives the operating temper- 
ature is such that the adhesive remains in its Initial glassy stage 
through the entire loading period and hence it is not necessary to con- 
sider the time-dependent behavior of the stress-strain relations in per 
forming the stress analysis of the joint. However, in certain appli- 
cations, the temperature and the load duration may be such that the 
rheological behavior of the adhesive may no longer be negligible. In 
this paper the adhesively bonded joint problem is considered by assum- 
ing that the adhesive is a linear viscoelastic material. 

2. Formulation of the Problem 

In formulating the adhesively bonded joint problems unless the 
thickness of the adherends is at least two orders of magnitude smaller 
than the length characterizing the bond region the generalized plane 
stress or the membrane assumption does not seem to be very realistic. 



( 


On the other hand In an adhesive joint between relatively thin adherend's, 
even if it were possible to formulate the problem by assuming the adher- 
ends as elastic continua, the numerical analysis Involve such severe 
convergence problems that the accuracy of the results may be highly 
questionable [4], In such problems the plate assumption in modeling 
the adherends appears to be a fairly good compromise. Thus, In this 
paper the problem will be formulated under the following primary assump- 
tions: (a) the adherends are treated as linear elastic plates and the 

transverse shear effects are taken Into account; and (b) the adhesive 
is considered as a viscoelastic solid in which the in-plate strain as 
well as out-of-plane strain and shear strain are assumed to be nonzero. 
The secondary assumptions under which the specific problem is formulated 
and solved simplify the analysis quite considerably but do not affect 
the character of the solution. These assumptions are: (a) the problem 

is one of plane strain, that is, the bonded joint is very "wide" and 
undergoes cylindrical bending; (b) the adherends have the same thick- 
ness and are made of the same material; and (c) the structure is a 
single lap joint. The elastic version of the problem neglecting the 
transverse shear effects in the adherends was considered in [2]. The 
solution of, again, the elastic problem for different adherends with a 
somewhat simpler adhesive model may be found in [S]'(*) ** 

The geometry of the problem under consideration is shown in Figure 
1(a). From the equilibrium of the plate elements for the adherends 1 
and 2 the following differential equations may be obtained: 
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(1 a-c) 


(2 a-c) 


where are respectively the membrane, transverse shear, 

and moment resultants, the index 1 = 1,2 referring to the adherends 1 


(*) Needless to say, the problem has been very widely studied. Some 
references to further analytical work and to finite element type 

solutions may be found in [3]. 
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and 2, h,, h«, and h are the thicknesses of the adherends and the adhe- 
sive as shown, and a(x,t) and T(x,t) are the interface normal and shear 
stresses. In modeling the adhesive it is assumed that the stress com- 
ponents a„(x,y,t) * a(x,t) and t (x,y,t) = T(x,t) in the adhesive layer 
y 

are independent of the y coordinate. 

Assuming cylindrical bending, = 0, ^ 2 ^ - 0* The stress resultant 
displacement relations may then be -expressed as 


3u, 36. 3v, 

i;r=ClNlx‘^= “l^lx’lT^^ 
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where 


1-Vj^ I2(l-Vj^) c 

, D. *» c k 3 » ^ ^i^j * ('“^2) 


E.h 
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(3 a-c) 
(A a-c) 
(5 a-c) 


E|» yj» Vj, (i=l,2) are the elastic constants, u.(x,t) and Vj(x,t), 

(i*»l,2) are x and y-components of the displacement vector and (1*1 >2) 

is the rotation of the normal to the midplane of the adherends. 

It may be seen that as stated the problem has lA unknown functions, 

name ly, d, t, ^i’ ^ix* ^Ix* ^ix* ^Ix* (*“^>2). Eguat i ons ( 1 — A) 
provide 12 relations. The remaining two relations necessary to complete 
the formulation of the problem are obtained from the continuity condi- 
tions for the displacements in the bond region. To do this the mechanics 
of the adhesive layer, specifically its constitutive relations need 
to be considered. 

Referring to Figure lb the average strains in the adhesive may be 
expressed as 

hi h2 

Txy “ ^“l “ T ^Ix " “2 " T ^2x^''^o ’ 


Sy = (v, - V2)/h^ , 

= 3Ui 36,^ 3 u2 

®x ~ ' 3x 2 3x 3x 2 3x 


(6 a-c) 
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Noting that all the remaining strain components in the adhesive are 
zero and defining 

e = (e^ + ey)/3, (7) 

the strain tensor for the adhesive may be decomposed as follows: 




1 

0 
O 

1 



0 0 0 


0 e 0 

0 0 e 

+ 

_ 0 0 -e 


Similarly, noting that 'o - a, t = t, the stress tensor for the 

y xy 

adhesive may be decomposed as 



where, the hydrostatic component of the stress tensor s is defined by: 

s = (a^ + 0 + a^)/3 (10) 

The constitutive equations of linear isotropic viscoelastic mater- 
ials may be expressed in terms of either hereditary integrals by using 
creep compliance or relaxation functions, or differential operators* [6-8]. 
In this paper the latter approach is adopted and it is assumed that 


P,(s.j) = Q,(e.j.), (i,j) = 1,2,3, 

(11) 

P2(s) = Q2(e) 

(12) 


(*) The two formulations are, of course, related through Laplace trans- 
forms. For example, the_creep compliance J(t) is the inverse Lap- 
lace transform of P(s)/sQ(s) where P and Q are the related differ- 
ential operators operating on a and e, respectively and s js the 
transform variable. 
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where s.j and e.j (i,j) = 1,2,3, are the deviatoric components of 
stress and strain tensors, respectively, as given by (8) and (9), s 
and e are defined by (10) and (7), and P^, P^, and are differential 

operators of the form S a,^(t) 3*^/3t'^, the coefficients a^ being generally 
functions of temperature. More explicitly, from (7)"(12) it may be 
seen that 


P,(2o^ - 0 - (Jj) - - Ey) , 

(13) 

P,(2o - cr^ - 0 ^) - a,(20y - e^) , 

(H) 

P,(2c7^ - "x - * V ■ 

. (15) 

P,(t) - ja,(Y^y) . 

(16) 

^'‘’x • «2<=x * • 

(17) 


3 3 

Since Z s.. = 0 and Z e.. =0, equations (13"15) are not linearly inde' 
I 1 1 T * * 

pendent. Equation (14) may be obtained by adding (13) and (15) and 
will, therefore, be ignored in the remainder of the analysis. 

Practical experience indicates that under a hydrostatic stress 
state most viscoelastic materials behave elastically. Hence, it may 
be assumed that 


P^ = 1, H 3K , (18) 

or 

+ a + = 3K(e^ + Cy) (19) 

where K is the bulk modulus of the adhesive. Eliminating and from 
(13), (15) and (19) and using (6 a-c) , the constitutive equations may 
now be written as 
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( 20 ) 

( 21 ) 
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Equations (20) and (21) with (1-4) provide the system of 14 relations 
necessary to solve for the unknown functions a, t, Uj, v., 0.^, N.^, 
and Mj^, (1=1,2). 


3. Example 


As an example we consider a single lap joint which consists of two 
identical adherends bonded through an adhesive layer which may be repre- 
sented by a three-parameter viscoelastic solid (Figure Ic). For the 
adherends we have 


h 


1 



h, 



D, = D, 


n ■_ 12(l-v^) 

" eP 


B, = B, 


5 

Z 


U h. 


( 22 ) 


For the adhesive, referring to Figure Ic it may be shown that 


^1 “ ' ®1 3t ’ ^1 “ •’o ^1 it 


(23 a,b) 


where 


^ 1^2 ^ 2^1 

*1 “ * ‘’o " " 1^2 


(24) 


For a nondecreasing strain under sustained load the following Inequality 
must be satisfied: 

b, > a^b^ . (.25) 

Generally, the coefficients ap b^ and b^ are functions of temperature, 
hence implicitly functions of time, if the temperature does not remain 
reasonably constant during the period of loading. In the example con- 
sidered, it is assumed that these coefficients are constant. 

I 

Through a relatively straightforward elimination, the goverryng. 
equations (1-4), (20) and (21) can be reduced to a pair of differential 
equations in the unknown functions o(x,t) and x(x,t) carrying out 

this elimination, using (22) and the operators defined fxy (.23) we obtain 



3^t . 3^t _ 

3x2 ^1 8x^3t 2h 


° [2C {h + h„)]T 

o ^ 


hn b^hD bihD 

2{^ [2C + -j- (h+ho)] ^ = “ - -Tjj— , 


3|#-^ 3a, ^ + [ f0(3K-b^) - j^(3K+2b,)] g 


3^0 


3^0 


( 26 ) 


+ Cf0(3Ka,-b,) - h^(3Ka, + 2b,)]^ 


+ D(3K+2b^)0 + IP D(3Ka,+2b,) |f = 0. 

o o 


(27) 


Assuming that no external transverse shear load is applied to the composite 
plate in -£ < X < £ and noting that r(x,t) is the average shear stress act' 
ing on the adhesive, referring to Figure la the equilibrium of transverse 
shear resultants gives 

Q,^(x,t) + Q2j^(x,t) + h^T = a^(t) , (-l<x<l) . (28) 

Equation (28) has been used in deriving (26). 

The differential equations (26) and (27) are uncoupled and may easily 
be solved by first reducing them to ordinary differential equations through 
the use of Laplace transforms defined by 

cx> 

F(s) = I f(t)e“®^ dt , (29) 

o" 

C+ i 00 

f(t) = 2I7 j F(s)e®^ ds (30) 

c-ioo 


where F(s) is called the Laplace transform of f(t) and the constant c is 
selected In such a way that all the singularities of F(s) lie to the left 
of the line of integration Re(s) = c. Assuming that the bonded joint is 
initially stress-free, the functions o(x,t) and T(x,t) are zero for t<0 
and from (26) and (27), we find 
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(31) 


d^F 


dx^ 


- a^F = e 




oj'^G = 0 


(32) 


where F(x,s) and G(x,s) are the Laplace transforms of T(x,t) and crCx,t), 
respectively and 


(Ac + hO(h+h^) ] (b^ + b.|S) 
Ah^d + a^s) 


(33) 


hDQ^(b^ + b|s) 
Ah^sCl+a^s) 


(3A) 


fh? “ (3K - b„) 

+ s[^ (3Ka^ + 2b^) - I D(3Ka, - b^)]} , 


(35) 


“ 3h; -f l ^ ^ iy ^ =(3Ka, + 2b,)] . (36) 

In the example it is assumed that the external loads are given by (see 
Figure la) 

No(t) = N^H(t), M,(t) = M,H(t), M2(t) = M2H(t), Q^(t) = Q^H(t) 

(37 a-d) 

where H(t) is the Heaviside function. For example, the nonhomogeneous term 
6 which appears in (31) and which is given by (3A) Is obtained by using 
(37d) . 

The general solution of (31) and (32) may be written as 

F(x,s) = sinh(ax) + A2 coshCax) - , (38) 

G(x,s) = sinh((j),x) + cosh((|),x) + A^ sinh(<(i2x) 

+ Ag cosh((|)2x) (39) 
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where 


<^1 = Cy ^ + . <^2 ' 

and the, unknown functions A^(s), ... , A^(s) are determined from the boun- 
dary conditions. 

The problem is solved under three separate loading conditions shown 
in Figure 2. 

(a) Membrane loading (Figure 2a). 

For this case the boundary conditions for plates 1 and 2 are given 
as follows: 

N,^(£,t) = 0, M^^(£,t) = 0, Q,^(£,t) = 0, (M a-c) 

N^^(.-£,t) = N^H(t), M^^(-£,t) = -N^ ^ H(t), Q^^C-£,t) = 0, 

(42 a-c) 

h+hf, 

N2^(£,t) = N^H(t), M2^(-f,t) = -j- H(t), 

(43 a-c) 

N2^(--e,t) = 0, M2j^(--e,t) = 0, Q2 x^“'^»^^ = ° » 


Considering the symmetry of the problem in geometry and materials, after 
some lengthy manipulations it can be shown that (41-44) are equivalent to 
the following conditions: 


x(x,t) 



dx = - N^H(t) , 


(45 a,b) 


a(x,t) = a(-x,t), J a(x,t) dx = 0 , (46 a,b) 

-I 


^(3K-b^)a(£,t) + |o(3Ka^-bp ^ o(£,t) - f^(3K+2bo)[ | a{l,t) 

+ f(h+h^)N^H(t)] - f^(3Ka^+2bp[|^ a(£,t) + Y(h+ho)N^6(t)3 
+ 3 a(£,t) + 3a^ alFFt = 0 . (47) 
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In this problem since = 0, 3=0 and substituting from (38) and 
(39) into (^5"^7) we obtain 


A,(s) = 0 , 

A 3 CS) = 0 , 

AjCs) = 0 , 


“ Nq 

" " irTTniTlair ’ 


A^Cs) 

AgCs) 


(h+hQ) Nqu‘*s i nh (.(ji2-^) 

4sT^aJ(s1 

(h+ho) Nqoj^^s i nh (({) ]£) 

lisT^TTs) 


m 


where 

Ag(s) = (|)2COsh((ti^£)sinh((j)2'^) - 4)jSinh(iJi^£)coshC(j)2-t) . (^9) 

(b) Bending (Figure 2b) 

For this problem the boundary conditions are 


N^^(£,t) = 0, M,^(£,t) = 0, ()^^(£,t) = 0 , (50 a-c) 


N,^(-£,t) = 0, M,^(-£,t) = M^H(t), Q,^(£,t) = 0 , (51 a-c) 


N2j^(-e,t) = 0, M2^(£,t) = M^H(t), “ ° » 


(52 a-c) 


N2x(--^.t) = 0, M2x(--e,t) = 0, ° ’ ^^3 a-c) 


Again, considering the symmetry of the problem conditions (50-53) may be 


shown to be equivalent to the following: 

T(x,t) = -T(-x,t) , (5^) 

't ^ ®1 ^ “ ■ •5K‘o^o[boH(t) + b,6(t)] . (55) 

a(x,t) = -a(-x,t) , (56) 

rl 

a(x,t)x dx = M H(t) , (57) 

j ^ 

-I 
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(3K-b^) ^a(£,t) + (3Ka^-b,) i^-^a(£,t) 

- ^^(3K+2b^) [| aU.t) + D MgH(t)] 

- i (3Kaj+2b,) + D 6(t)] 

+ 3 aU.t) + 3a, cr(£,t) = 0. (59) 

In this problem, too, 6*0, and substituting from (38) and C39) into 
(5^*59) we obtain 


hDMo(bo+b,s) 

ijh^as (1+a , s ) cos h (a£) * ^2^®^ ® ’ 


A 2 (s) 


cosh((}> 2 -C) 
2 s (^2 A[j(s) 


A 5 (s) * 


coshUiZ) 

2 s 4 >,Ajj(s) 


Ai,(s) = 0 , 

Ag(s) * 0 


A|j(s) * c{) 2 sinh((i),^)cosh(((» 2 -£) " <|), cosh ((}>,£) s i nh C«j> 2 ’^) • (^O) 


(c) Transverse Shear (Figure 2c) 

For the loading given in Figure 2c the boundary conditions may be 
expressed as follows: 

N,^(-e,t) * 0, M,^(£,t) = 0, Q,j^(-£,t) = 0 , (.61 a-c) 

N,^(-£,t) = 0,M,^(-£,t) = -Q/H(t), Q,^(-£,t) - O^H(t) , C62 a-c) 
N2x(-e,t) * 0, M2x(-e,t) = Qo^H(t), , (63 a-c) 

N2x(--^»t) * 0, M2x(“4t) = 0, d2x(“'^*^^ “ °* 

These conditions are equivalent to 

(I , 

t(x,t) = T(-x,t), I x(x,t) dx * 0, C65 a,b) 

-Z 
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C66 a,b) 


a(x,t) = aC-x,t), j c(x,t)dx = - Q^HCt), 

-I 

' ^(3K-bQ)a(£,t) + i^(3Ka^-b,)^ a(£,t) - i^(3K+2b^)[| aU,t) 

+ DQ^£H(t)] - -J-^(3Ka,+2b^)[|^ a(£,t) + DQ^£5(t)] 

+ 3 ^ aU,t) + 3 a, aU,t) = 0. (6?) 

In this case S is given by (3^) and, the functions (s) , . . . ,A^Cs) are 
obtained as follows: 

■A,{s) = 0, k^(s) - „ sffhtat) ’ 

AjCs) = 0, A^(s) . '^[“^'l>lcosh(t?a -^^>£slnh(^;,£)] _ 

A (s) = 0, A (s) = - ~ sinh((|)i£)1 

5 ’6 ^ 2 s <j) 1 (s) ’ 

A^(s) = (|) 2 Cosh(({>^£)sinh((J) 2 '£) “ sinh((j)^£)coshC'|) 2 '^) • (-68) 


A. Solution and Results 

After determining the functions Aj (s) , (1=1,..., 6) the unknown func- 
tions t(x,t) and a(x,t) may be obtained by substituting from (38) and 
(39) into the inversion integral (30). In each case the constant c giving 
the line of integration is determined by analyzing the singular behavior 
of the functions F(x,s) and G(x,s) in the complex s plane. Because of 
the existence of a number of branch points in the complex plane the exact 
inversion of F and G becomes very complicated and, in light of the fact 
that the inversion integrals can be evaluated in a straightforward manner 
numerically, does not seem to be worth the effort. Thus, making the 
following change in variable 

s = c+ iy, -®<y<co ( 69 ) 
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the functions t and a may be expressed as 

00 

T(x,t) F(x,c+iy)e^^‘^‘*’'^^dy , 

' —00 

o(x,t) G(x,c+iy)e^^‘^'*’‘^^dy , 

-00 

It can be shown that the imaginary parts of the integrands in (70) and (71) 
are odd functions in y and therefore the integrals give real results. 

Examining the functions F and 6 in- the complex plane it is found 
that s=0 is a simple pole and all the remaining singularities lie in the 
left hand plane. Hence c is a positive constant. To evaluate the inte- 
grals in (70) and (71) first they are expressed in (0,<») as follows: 


C70) 

C71) 


*00 *oo 

f(y)dy = [f(y) + f("y)]dy. (72) 

— oo 0 

Even though there are routine techniques for evaluating infinite integrals, 
it is generally a good practice to obtain the asymptotic behavior of the 
integrands for large values of the argument before selecting a particular 
technique. In the problem under consideration the integrands do not decay 
exponentially. Consequently, the numerical integration requires more 
care. One way to insure that no significant accuracy is lost due to the 
slow decay of integrands is to evaluate the integral in closed form for 
large values of the argument. For example. In the lap joint under mem- 
brane loading Nq (Figure 2a), after analyzing the asymptotic behavior of 
the function F, the shear stress t may be expressed as 


t(x,t) = 

.!l£ 

2tt 



r-f .. ^ t(c+iy) 
[F(x,c+iy)e 

cosh(x/m/a^) 
® sinh(£»/rn7a,T 


+ F(x,c-iy)e^^*^ ''^^Idy 


1 


sin(ty) 

y 


dy 


where 

m = b^[AC + h(h+hQ)D]/(i»hQ) 


(73) 

(74) 
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and A Is a "large" number. The second integral Is known in closed form 
and the first is evaluated numerically. The proper selection of A 
requires some trial calculations. In this problem A selected in 20 to 
30 .range gives good results. It may also be pointed out that the numer- 
ical calculations show the results to be insensitive to the choice of 
the constant c. 

in the numerical example considered it is assumed that the adherends 
are aluminum alloy plates with the following elastic constants and dimen- 
sions (Figure 2) 

E = 10^ psi = 6.895 X 10^° N/m^ , v = 0.3 
h = 0.09 in = 0.229 x lo‘^m , Z = 0.5 in = 1 .2? x lo'^m. 


In the three parameter viscoelastic solid adopted for the adhesive the 
coefficients which appear in the operators P| and Csee eqs. 23a, b) 
are related to the constants shown in Figure Ic by (24). To relate these 
constants to somewhat more conventional material properties consider the 
response of the model given in Figure Ic to an input t= XoUCt) which is 
found to be 


1 


(0 =^Ct.(l-e"^''^o) + a,e"^/^o], t. 


'1 



C75) 


where t is called the retardation time. Now defining 
o 


o o 

^o ^ ^ 

from ( 7 ^ it is seen that 

bi ^1^9 

2^0 = = '^1 = bo = 1^2 


C76) 


(77) 


Thus, the moduli u and u and the retardation time t may be selected 
as the three parameters representing the viscoelastic solid. 

For the particular epoxy used as the adhesive the properties at 
t=0 are assumed to be 
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-k 

h = 0.00^ in. = 1.016 x 10 
o 

E = 5.797 X 10^ psi = 39.968 x 10® N/m^ 

o 

= 2.225 X 10^ psi = 15.3^1 x 10® N/ra^ 


The bulk modulus K is assuned to be constant and may, therefore, be cal- 


culated in terms of and the shear modulus as 

^0^0 


1/ — ^ ^ 

^ 3(3y -E J 


(78) 


o o 


In the example it is also assumed that 


y^ = y^3 » ^ hrs. 

If it is assumed that the adhesive layer is linearly elastic having; 
the constants E and v , with the adhesive model used in this paper the 
solution may be obtained in a straightforward manner. For example, in the 
case of membrane loading described by (A1-A4) the adhesive stresses are 
found to be 

cosh (a^x) 

. (79) 


T,(x) =- 


o e ^ e 

2 sinh (a £) 


“e^ “ 4(l4 ~ )h ~ + hD (h + h^)] , 


a' o 


( 80 ) 


a^(x) = B^cosh(m^x) + BgCosh(m 2 x) 


( 81 ) 


B, = - e.^^N (h+h )sinh(m,£)/(Am,A) , 
loo ' ' 


I — — c. , »'l V I I J I I II* ^ 

H loo Z ^ 

B, = (h+h )sinh(m,£)/(Am,A) , 

o 1 O O 1 I 

A = m2COsh(m|£)sinh(m2-£) - m^sinh(m^£) cosh(m2-&) . 


*^1 h„0-vV) ’ 


0-v.)£. 


2 _ 'a' a r 1 


hOv. 


" (l-v^-2v^^) Bh ■ A(l-v ) 

3 3 U a 


] . 


( 82 ) 




J? 
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On the other hand, in the case^viscoelastic adhesive the elastic 
response for t=0+ and t=« may also be determined by using the limit theorems 
for the inversion of Laplace transforms. For example, again for the case 
of -membrane loading, from (38) and (A8) the shear stress in the adhesive 
may be obtained as 


N a cosh (a x) 

- - 2 sinh - 

o 1 

N a cosh (a x) 

2sinh(aJ) ’ 

b 

®oo^ ~ 2j(^ ^'^(^■*'^o)l • 


C83) 

(81») 

(85) 

C86) 


Note that at t=+0 p =E /2(l+v,) = and = E^, and from C77) , (.80), 
and (8A) it follows that = a^. Hence, the initial response given by 
(83) is the expected elastic solution given by (79)* Similarly, at t=» 
p =p , and (77), (80) and (86) shows that = a., and hence t(x,~) = 

Tg(x). Also, it can be shown that a(x,“) corresponds to the elastic solu- 
tion obtained by using p =p and the bulk modulus of the adhesive which 
is assumed to be a time- independent constant. 

For the three types of loading shown in Figure 2 the calculated results 
for r(x,t) and a(x,t) are shown in Tables 1-6. To visualize the variation 
of the stresses in time and along bond region some sample results are also 
given in Figures 3-5. Figure 3 and k show the distribution of shear and 
tensile stresses in the bond region in a single lap joint under membrane 
loading for some fixed values of time. As expected, there is a certain 
redistribution of stresses with increasing time. This may also be seen in 
Figure 5 where the variation of the maximum values of t and a is given. 

From Figures 3-5 and Tables 1-6 it may be observed that the peak values of 
the tensile stress a in the adhesive are not only higher than the correspond- 
ing shear values but also decay slower. The values x and a given in Tables 
1 and 2 for t=0 and t=“ are obtained from the elastic solutions (.79) and 
(81) by using the bulk modulus K which is assumed to be independent of time 
and the corresponding p^ and p^. 
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Table 1. Variation of x(x,t)/(No/£) for 

the case of tension (t in hours) 


x/Z ^ 

0 

0.01 

0.1 

0.5 

1 .0 

2.0 

4.0 

00 

0 


HIHSI 

BWHIM 




BSSSS^EB 

-0.012 

0. 1 

-h 

-3.67x10 

BBiBHI 

BiHeei 





-0.015 

0.2 

-If 

-9.81x10 

■■BUM 

-1.16x10"^ 



ffffiWHHI 

-0.010 

-0.023 

0.3 

-2.89xl0“^ 

rnSSSmSm 

BWfflfBI 



-0.013 

-0.022 

-0.041 

0.^ 




-0.014 

-0.020 

-0.030 

-0.048 

-0.075 

■0.5 


-0.025 

-0.027 

-0.038 

-0.050 

-0.071 

-0.101 

-0.139 

0.6 

-0.072 

-0.073 

-0.078 

-0.100 

-0.125 

-0.164 

-0.213 

-0.258 

0.7 

-0.212 

-0.213 

-0.224 

-0 . 266 

-0.309 

-0.373 

-0.441 

-0.481 

0.8 

-0.62^f 

-0.624 

-0.641 

-0.702 

-0.759 

-0.833 

-0.891 

-0.896 

0.9 

-1 .83^1 

-1.828 

-1.838 

-1.843 

-1.839 

-1.816 

-0.757 

-1 .670 

1.0 

-5.391 

-5.351 

-5.265 

-4.812 

-4.382 

-3.838. 

-3.366 

-3.112 














































Table 2. Variation of o(x,t)/(No/-£) for 
the case of tension (t In hours) 



1.01x10 


1.18x10 


0.1 l-^». 61x10 33x10 


0.3 -2.i»5xlO~-'|-2.1»4xlO 


0.1 0.5 1.0 


2.05x10"^ 5.80x10"^ 9.57x10"^ 


-3.01x10“^ 2.77x10"^ 8.59x10“^ 



W 
















































































Table 3» Variation of x(x,t)/(MQ/£^) for 
the case of bending Ct in hours) 
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Table 4. Variation of a(x,t)/(MQ/£^) for 
the case of bending (t in hours) 
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Table 5. Variation of TCx,t)/C0o/^) 

the case of shearing Ct in hours) 
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Table 6. Variation of aCx,t)/(Q^/£) for 

the case of shearing (t in hours) 
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Figure 5. Variation of maximum shear and normal stresses in the 
adhesive as functions of time. 
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